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Abstract. In the present work we approximate an ultrarelativistic jet by a homo- 
geneous beam of null matter with finite width. Then, we study the influence of this 
beam over the spacetime metric in the framework of higher-derivative gravity. We flnd 
an exact shock wave solution of the quadratic gravity field equations and compare it 
with the solution to Einstein's gravity. We show that the effect of higher-curvature 
gravity becomes negligible at large distances from the beam axis. We also observe 
that only the Ricci-squared term contribute to modify the Einstein's gravity predic- 
tion. Furthermore, we note that this higher-curvature term contribute to regularize 
the discontinuities associated to the solution to Einstein's general relativity. 
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1. Introduction 

The relativistic and ultrarelativistic jets arise as important components in the structure 
and dynamics involved in several astrophysical scenarios. Among the different types 
of astrophysical jets, the most energetic ones are potential candidates to give rise to 
emission of gravitational waves. For example, highly relativistic jets should be associated 
with some sources of gamma ray bursts (GRBs). In the last few years, observations 
that indicate the presence of jets in GRBs have been reported by many authors pp. 
The expansion of the jet in the burst leads to very high Lorentz factors which is at 
least of hundreds, but higher values such as F ~ 10^ — 10^ are taken into account in 
modelling the ORB sources [SJ. In the present paper, we consider a model in which the 
gamma Lorentz factor takes the limit T ^ oo. This assumption can be used to solve 
exactly the gravitational field equations. Then, we can take this analytic solution to 
represent the gravitational wave generated in a strong burst. To study the impact of 
an ultrarelativistic jet over the spacetime metric, we start from the extreme situation 
where the velocity of the particles in the beam is assumed to be equal to the velocity 
of light. The jet is then represented by a beam of null particles. For jets which start 
with a small opening angle Oq < 10"^ — 10~^ PP, we assume that the width of the beam 
remains constant during the first stage of the jet expansion. Then, we calculate the 
effect of this jet over the spacetime metric in a flat background. We take advantage of 
the use of exact methods to provide solutions to higher-derivative gravity equations 
and compute the solution to quadratic gravity field equations. 

Quadratic gravity is an example of the higher order theories of gravity which are 
generally covariant extensions of the general relativity. These theories are constructed 
by the addition of terms nonlinear in the curvature to the Einstein-Hilbert action. 
The coupling parameters of the quadratic curvature terms in the Lagrangian must 
be determined by experiments. Observational constraints on these parameters can be 
traced out from sub-millimeter tests of the inverse square law |3] or by the bending of 
light by a gravitational field 0. Are there a way to obtain an observational constraint 
on quadratic gravity coupling parameters from gravitational radiation? The application 
of the traditional methods of ordinary gravitational wave theories, which deals with 
the linearized version of the gravitational theories, are until the present date, unable to 
clarify the role played by the quadratic curvature terms in the gravitational radiation [H]. 
In a recent paper [Zj we deal with the linearized version of quadratic gravity treating 
the higher-curvature terms perturbatively. We conclude that the effects of quadratic 
curvature terms are so small that cannot be measured by the current interferometric or 
mass resonant detectors \. 

In the present paper, we obtain a solution which can be able to explain the effect 

I The amplitude of a general linearized oscillatory quadratic gravity wave at a given distance from the 
source differ from the Einstein's linearized wave amplitude by a frequency dependent function. For a 
wave with frequency of 100 Hz and considering an upper bound of lO^^'cm^ for \f3\ this difference 
becomes < 10~^^. This difference cannot be measured by the current gravitational wave detectors. 
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of higher-curvature invariants on gravitational radiation by deahng with exact solutions 
to gravitational field equations. We obtain exact shock wave solutions. A gravitational 
shock wave is truly a discontinuity of the spacetime metric which propagates with the 
velocity of light in a given direction jjSj. The notion of wave amplitude as a finite 
quantity which varies with the distance from the source cannot be applied here. The 
gravitational shock wave must be taken instead as a propagating singularity. According 
to [H] a modification of the gravitational wave detectors would be necessary to observe 
gravitational shock waves. We do not intend to discuss the detection of gravitational 
shock waves in the present work. The tidal field of a gravitational shock wave is an issue 
which we are dealing with and will be presented in a forthcoming paper. 

The plane of the present paper is the following. In section |21 we write the field 
equations for the shock wave metric in the frameworks of quadratic gravity and Einstein's 
general relativity. This is carried out by assuming that the solution can be described 
by a pp-wave metric jTUj . In section El we consider a homogeneous and finite beam of 
null particles and obtain the gravitational shock wave generated by it. We also compare 
the solutions obtained in both theories. In section |3] we summarize the principal results 
obtained and make some comments concerning to the application of these results in 
astrophysics. 

2. The field equations for a plane gravitational shock wave 

We start with an impulsive pp-wave metric given by the following line element 

ds"^ = -dudv + H{x\ u)du^ + {dx')^, (1) 

where 

H{x\u) = f{x')5{u), u = t-z, v = t + z, (2) 

and (i=l,2) denote the Cartesian coordinates in the plane perpendicular to the 
wave propagation, which we call transverse coordinates. The metric represents a 
gravitational shock wave propagating in z direction. This is a variation of the Peres's 
wave metric [TTj and was used by many authors to study wave-like exact solutions to 
Einstein gravitational equations |H1 The coefficient of ^-distribution is the wave 
profile function. This function must be determined by the field equations and will 
depend on the characteristics of the source and on the underlying theoretical model. 
The nonvanishing Christoffel symbols for the metric are given by 

Ku = -fs, Ku = -ldj6, ^^ = = -dj6, (3) 

where di denotes partial derivatives with respect to the coordinate x\ and the dot 
denotes the derivative with respect to the u coordinate. The only nonvanishing 
components of the Riemann tensor, apart from the ones obtained by symmetry 
properties, are given by 

Riuju = ~didjH{x\ u). (4) 
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The only nonvanishing components of the Ricci tensor for the metric (Q) are 

Ruu = -lvlH{x\u), (5) 

where V^^ denotes the Laplacian operator in the transverse space {x*}. The curvature 
scalar R and the elementary quadratic invariants vanish identically 



R = 0, RapR''^ = Rap^sR'"^^^ = 0- (6) 
The quadratic gravity theory in four dimensions can be derived from the action 

S = 77^ / (fx^{R + aR^ + ISR^.R^" + 167rG/:„}. (7) 
IdttG J 

Since, in four spacetime dimensions, the Gauss-Bonnet invariant can be used to eliminate 
the Ra/B'ysR"''^'^^ term. In the above action, stands for the presence of matter fields. 
The units are such that h = c = 1. We remark that the gravitational actions built 
from Lagrangian densities which are arbitrary functions of R are conformally equivalent 
to general relativity interacting with scalar fields. However, this equivalence cannot be 
traced for the theory derived from the action ((7j) due to the presence of the R^uR^'^ 
invariant |13j . 

Requiring the action S to be stationary leads to the following field equations 

Rfiu - ^gtMuR + aH^i. + PlfMu = SttGT^^, (8) 

where H^,^ and If^u are the operators associated respectively with the R^ and R^yR^" 
invariants and are given by 

H^, = -2R,^y + 2g^yUR - ^g^.R^ + 2RR^y, (9) 

and 

^tj,u = —'^R'^^-va + '^R^lV + -jdnv'^R + '^R'^fj.Rav — ■^QfiuRa/sR"'^ , (10) 

where □ refers to the curved space d'Alambert operator. The theory defined by the 
action (I7|) leads to the following non-relativistic gravitational potential |14j : 

— + \- , (11) 

r 6 r 3 r J 

where 

This potential reflects the three modes of the linearized theory, one with a vanishing 
mass which gives the Newtonian force and two massive modes which create Yukawa-type 
interactions. To obtain an acceptable Newtonian limit we restrict mo and mi to real 
values, leading to the, so called, no-tachyon constraints 

3a + /3>0, -/3>0. (13) 
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Experiments which tests the inverse square law at sub-milhmeter distances indicate 
an upper hmit to the absolute value of the a or j3 parameters of 10~^cm^ A direct 
observational constraint on the f3 parameter can be obtained from the bending of light 
by the Sun's gravitational field. According to a semiclassical computation of the bending 
of light in quadratic gravity 0, the higher-curvature coupling parameter (3 must satisfy 
the same constraint, namely |/3| < 10~^cm^. This bound improves by several orders of 
magnitude the early accepted value for the upper limit on [T3] . 

Now, let us write the field equations for the metric given in (^. By virtue of the 
identities © the operator H^j, vanishes, ensuring that the square curvature scalar does 
not contribute to the impulsive wave solution. The field equations (jHl), reduces to the 
following linear fourth-order two dimensional partial differential equation: 

- + ^l]H{x\ u) = 87rGT„„, (14) 
where = V5_V5_. In the Einstein gravity, a = P = 0, the equation (jSI) reduces to 

- I'^lHEix', u) = 87iGTuu{x\ u). (15) 
Taking into account the above equation we can rewrite (jl4|) as 

[/?Vi + l]H{x\ u) = He{x\ u) + ah{x\ u), (16) 

where h{x^, u) is a harmonic function of the transverse coordinates and, a an arbitrary 
constant. This result means that we can add to the final solution a harmonic function 
of the transverse coordinates. The particular dependence of h on is dictated by the 
symmetries of the problem and the dependence on the u coordinate is determined by 
the source term. By making a decomposition of H{x\ u) as 

H{x\ u) = Hq{x\ u) + He{x\ u) + ah{x\ m), (17) 

where the index Q refers to the purely quadratic part of the solution, we obtain the 
following second order partial differential equation for Hq{x'^,u): 

Hq{x\u) = 16nGT^u. (18) 

Thus, the problem of solve the fourth order equation (fT^ was reduced to the problem 
of solve the second order equations (fT3j) and (fTH|). 



V + 



P. 



3. The gravitational shock wave associated with a homogeneous beam of 
null particles 

In this section, we compute the solution of equations (fT3jl and (fTH|) when the source is 
given by a distribution of null particles moving along the same direction. The relevant 
component of the corresponding energy momentum tensor in polar coordinates is given 
by iSi: 

Tuu = \0ir,9)5iu), (19) 
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where g{r, 6) is the energy density of the distribution of null particles and A is a 
dimensionless constant which varies within the range 1 < A < 2 and comes from the 
equation of state 

P = {\-1)q, (20) 

which describes an ultrarelativistic perfect fluid [THj . 
The line element ^ can be written as 

ds^ = -dudv + [/E(r, 6) + /Q(r, e)]6{u)du^ + dr^ + r^dO^, (21) 

where the functions /e and /q must satisfy the following equations: 

VVE(r,^) = -167rGAf?(r,e) (22) 



and 



Ur,9) = 167rG\gir,9). (23) 



We assume that the energy density is a constant qq within a certain radius 
< r < Rq and zero outside. Thus, the source represents a cylindrical beam of null 
particles with width Rq. This beam is a simple generalization of a single null particle. 
We know that the wave profile of a gravitational shock wave generated by a single null 
particle in quadratic gravity is given by 



/o(r) = 8Gp 



ln{-]+Ko 



(24) 



where p is momentum of the particle, Kq a modified Bessel function and tq some 
integration constant jT2|. This is a continuos function of r which is regular at the 
origin an diverges logarithmically at r ^ oo. 

To obtain a solution which generalizes the gravitational shock wave generated by 
a single null particle we proceed by solving the equations and (PHj) choosing the 
integration constants in such a way that the wave profile /(r) = fE^r) + fq^r) fulfill 
the same boundary and regularity conditions that the single null particle shock wave 
profile ()24j) . Therefore, /(r) must be i) a continuous function of r ( even at r = Rq, 
where the source have a discontinuity ), ii) regular at the origin, namely /(O) < oo, and 
Hi) logarithmically divergent at r oo. 

Let us solve the equation ()22|). The cylindrical symmetry implies that /e will 
depend only on the r coordinate. It is easy to see that the solution to (j22|l can be 
written as 

f -AnGXgor' + Ci ln(r) + G^ for r < Rq 

[ Di ln(r) + D2 for r > Rq. 

where Ci, G2, Di and D2 are integration constants. 

By intergrating both sides of (f^ . with /e(?") given by (j^ . over a circle of radius 
r < Rq and applying the divergence theorem to the left-hand side term we get Ci = 0. 
Now, integrating over a circle of radius r > Rq and applying the divergence theorem 
to the left-hand side term we obtain Di = —SnGXgQRl. The remaining integration 
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constants C2 and D2, are determined by the boundary and regularity conditions. As a 
regularity condition we impose the continuity of /e('") at r = Rq. This condition imply 
that 



D2 = IGnGXgoRl Q ln(i?o) - ^ + C^. 



The constant C2, which gives the value of /e at r = 0, is the only arbitrary constant that 
will be determined by the boundary conditions. The regularity at the origin is satisfied 
by given any finite value to C2, Then, for simplicity and without loose of generality we 
choose C2 = 0. Thus, the solution of equation for any r > is given by 



SuGXqqR^ 



In 



+ 



for 
for 



r<Ro 
r > Ro. 



(26) 



Let us turn out to the purely quadratic gravity part of the wave profile, /Q(r) 
Integration of the equation leads to the following result [TH] : 

1,2 



A,Ko ('^) + A2I0 Q - l67rGXgob-' 



for 



r<Ro 



(27) 



for r > Rq, 
'—(3 and Ai, A2, Bi and B2 are 



where ly and Ky are modified Bessel functions, h e 
integration constants. 

To obtain a solution /(r) regular at the origin and logarithmically divergent at 
r ^ 00, we set S2 = and = in (j^Zj). As we want that /(r) be a continuous 
function of r, /Q(r) must be continuous at r = Rq. This condition is fulfilled by taking 



A2 = WTiGXgobRoKi 



and 



5i = - WnGXgobRoIi 



Rn 



Thus, fqir) takes the following form 



fQir) 



IGTrGXgo 



RobKi 



Rq 
T 



for T < Rq 



(28) 



-IQ^GXgMh Ko Q for 

The solution to equation (fT^ with the source term (fT^ is given by 

H{r,u) = [f{r) + ah{r)]6{u), 



r > R 



0, 



(29) 



where /(r) = /£(?") + /q('") is the wave profile and h{r) is a harmonic function. The 
cylindrical symmetry implies that h{r) = ln(r). Note that h{r) diverges at r = and 
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Figure 1. A comparison between ttie wave profiles in Einstein's gravity (dashed line) 
and for quadratic gravity (solid line) . To clarify the difference of the predictions given 
by Einstein and quadratic gravity, we set AnGXgo = 1, Rq = 1 and plot the /(r) curve 
for \/—(3 = 0.3. We note a non zero residual value at r = for the quadratic theory 
and, we als note that the two curves becomes indistinguishable after a few Rq distance. 



r = oo. If we choose a = —SttGXqoRI the wave profile approaches to a constant value of 
a/2 as r — > oo. However, the solution will diverge at r = unless that we put Ai = —a. 
But, this choice violate the previous choice, Ai = 0, and implies in the lack of continuity 
of fQ{r) at r = Rq. Therefore, to fulfill the boundary and regularity conditions chosen 
for /(r) we must set a = in ()29j) . The final result diverges logarithmically at r oo 
and resembles the wave profile for a gravitational shock wave generated by a single null 
particle [12] 

The figure HJ shows the r dependence of the wave profiles fE^r) (dashed line) and 



This quantity results from the Ricci-squared term in the gravitational Lagrangian. It 
carries the information of the beam width Rq and tends to zero as Ro — »• 0. We also 
observe that the wave profile obtained for quadratic gravity agrees with the prediction 
of Einstein's gravity at large distances from the beam axis. 

A comparison between the first and second derivatives of the wave profiles obtained 
in the two theories shows that the higher curvature terms of quadratic gravity 
contributes to smooth out the r dependence of these functions at r = Rq. This 
"smoothing" property of quadratic gravity is related to the fact that the field equations 



f{r) = /E(r) + /Q(r) (solid line). 

The non zero value of /(r) at r = is given by 




(30) 
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Figure 2. In this figure we plot f'{r) (solid line) and f^{r) (dashed line) as function 
of the polar coordinate r. The units are the same used to plot figure ^ but here we 
set /3 = 0.07 to compute /'(r). We note the softness of the quadratic gravity result 
at r = i?o in comparison to the Einstein's gravity one. 
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Figure 3. Here are depicted the f"{r) (solid hne) and f^{r) (dashed hue) as function 
of r using the units of figured with /? = 0.07. We note that the quadratic gravity 
curve is smooth at r = _Ro while the Einstein's gravity curve shows a discontinuity at 
this point. 

are of fourth order in quadratic gravity, whereas they are of second order in Einstein's 
Gravity. 

In figure |21 we plot the first derivative with respect to r of the wave profiles in 
both theories. The dashed line represents f^ii^) and the solid line represents f'{r) (the 
prime denotes the derivative with respect to r). Figure El shows the r dependence of 
the second derivatives of the wave profiles. We note the continuity and softness of f"{r) 
at the point in which f^{r) has a discontinuity. Therefore, we can conclude that the 
quadratic curvature component associated to the Ricci-squared term in the gravitational 
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action smooths out the functional dependence of the derivatives of the wave profile. 
Furthermore, this higher-curvature term contribute to remove the discontinuity of the 
second derivative of the Einstein's gravity wave profile at r = Rq. 

We must bear in mind that a gravitational shock wave cannot be treated in the 
same way as a linearized gravitational wave. The effect of a gravitational shock wave 
on test particles is a subject that we will treat in an incoming paper. Nevertheless, we 
can remark that the test particles would be be sensitive to the analytic proprieties of 
the wave profile function JJj- Obviously, these properties remains the same whatever 
the particular values assumed for the (3 parameter. Thus, even at the limit of 
the quadratic gravity could in principle be distinguished from Einstein's gravity. 

4. Final remarks 

Let us now summarize the main results obtained in our investigation and give a 
suggestion for application of these results in astrophysics. 

The first remarkable result is that the quadratic invariant does not contribute to 
the wave solution. The only contribution, which comes from quadratic part of the theory, 
comes from the Ricci-squared term in the quadratic gravity action. The effect of this 
term over the gravitational wave emission is to regularize the discontinuity associated 
to the Einstein's gravity solution. We recall that at large distances from the beam axis 
the predictions derived in both theories coincide. 

The gravitational shock wave is very different as compared to an ordinary 
gravitational wave ^Sj. The geodesies and the geodesic deviations in spacetimes of 
gravitational shock waves have been studied in [I7j. The authors have showed that 
the behavior of small geodesic deviations will be dependent on the combinations of the 
wave profile and theirs derivatives in the transverse space with products and powers of 
the Dirac 6 distribution and the "kink" function, defined by u6{u), where 6{u) is the 
step function. Therefore, the effect of a gravitational shock wave on test particles will 
be dependent on these results. We consider that the issue of observation of gravitational 
shock waves deserves a separated investigation and we intend to do so in another paper, 
now in preparation, to appear elsewhere. 

An important question remains to be discussed. Can the model studied in present 
the work be addressed to the study of astrophysical sources? We know that some 
GRBs are generated in the expansion of ultrarelativistic jets QUI]- Then a cylindrical 
beam of null particles can be used as a frist approximation to an ultrarelativistic jet. 
In this approximation an exact shock wave solution was obtained. We compute the 
wave forms in the framework of Einstein's and quadratic gravity. If a gravitational 
shock wave takes place in a GRB scenario it may have important consequences to the 
physical process at the source. According to [THj a gravitational shock wave would 
be relevant to the high energy scattering of quantum particles. Therefore it can be 
espected that the gravitational shock wave contributes to the dynamics of a GRB 
source and to the composition of the ejected material. Furthermore, once a gravitational 
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shock wave related to a GRB can be detected, the results derived in the present work 
could contribute to the understanding of the ultrarelativistic jets features from their 
gravitational shock wave profiles. 



Acknowledgments 

We thank Dr. J. A. de Freitas Pacheco for stimulating discussions. We also thank 
the Brazihan agency FAPESP (grants 00/10374-5, 97/13720-7, 97/06024-4) and CNPq 
(grants 300619/92-8, 380.503/02-6) for the financial support. 



References 



[1] Sari R, Piran T and Halpern J P 1999 Ap. J. L17 519; Sari R 2000 Beaming and jets in gamma 
ray bursts, astro-ph/0002165, Piran T 2000 Phys. Kept. 333, 529-553; Meszaros P 1999 Origin 
of gamma ray bursters, astro-ph/99 12546 

[2] Usov V 1994 Mont. Not. R. Ast. Soc. 267, 1035; Usov V and Smolsky M 1997 Wide ultrarelativis- 
tic plasma beam-magnetic barrier collision and astrophysical applications, astro-pli/9704152 
Rees M J and Meszaros P 1992 Mont. Not. R. Astr. Soc. 258, 41. 

[3] Buchdahl H A 1983 J. Phys. A: Math. Gen. 16, 1441. 

[4] Long J C, Chan H W and Price J C 1999 Nud. Phys. B 539, 23-34, |h ep-ph/9805217| Hoyle C D, 
Schimidt U, Heckel B R, Adelberger E G, Gundlach J H, Kapner D J and Swanson H E 2001 
Phys. Rev. Lett. 86, 1418, hep-ph/0011014" 

[5] Accioly A and Bias H 2001 Phys. Rev. D 64 067701, gr- qc/0 1070*031 

[6] Will C M 1993 Theory and Experiment in Gravitational Physics (New York: Cambridge 
University Press). 

[7] de Rey Neto E C, Araujo JCN and Aguiar O D, Class. Quantum Grav., submitted. 

[8] Barrabes C and Israel W 1991 Phys. Rev. D 43 N. 4, 1129. 

[9] Barrabes C and Hogan P A 2001 Int. J. Mod. Phys. D 10, 711-722. 
[10] Kramer D, Stcphani H, Hcrlt E and MacCallum M 1980 Exact Solutions of Einstein's Field 

Equations (London: Cambridge University Press) p 233. 
[11] Peres A 1959, Phys. Rev. Lett. 3 N. 12, 571. 
[12] Campanelli M, Lousto C O 1996 Phys. Rev. D 54, 3854. 

[13] Gottlober S, Schimidt H J and Starobinsky A A 1990 Class. Quantum Grav. 7, 893-900. 
[14] Accioly A, Azeredo A, Mukai H and de Rey Neto E C 2000 Prog. Theo. Phys. 104, 103. 
[15] Stelle K S 1978 Gen. Rel. and Grav. 9, 353-371; Donoghe J F 1994 Phys. Rev. Lett. 72, 2996. 
[16] Neilsen D W, Choptuik M W 2000 Class. Quantum Grav. 17, 733. 

[17] Kunzinger M and Steinbauer R 1999 J. Math. Phys. 40, 1479; Steinbauer R 1998 J. Math. Phys. 
39, 2201. 

[18] Gradshteyn I S and Ryzhik I M 1965 Table of Integrals, Series and Products^ (New York: 
Academic Press). 

[19] Garriga J and Verdaguer E 1991 Phys. Rev. D 43, 391; 't Hooft G 1987 Phys. Lett. B 198, 61; 



Lousto C O Sanchez 1994 Scattering processes at the Planck scale, gr-qc/9410041 



